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£\j ^ Abstract. We review the Reidemeister torsion, Ray-Singer's analytic torsion and the Cheeger- 

Miiller theorem. We describe the analytic torsion of the de Rham complex twisted by a flux 
form introduced by the current authors and recall its properties. We define a new twisted 
analytic torsion for the complex of invariant differential forms on the total space of a principal 
circle bundle twisted by an invariant flux form. We show that when the dimension is even, 
such a torsion is invariant under certain deformation of the metric and the flux form. Under 
T-duality which exchanges the topology of the bundle and the flux form and the radius of the 
circular fiber with its inverse, the twisted torsions of invariant forms are inverse to each other 
for any dimension. 
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00 . Introduction 

Cn| ' Reidemeister torsion (or i?-torsion) was introduced by Reidemeister |26j for 3-manifolds. It 

. was generalized to higher odd dimensions by Franz [9] and de Rham [7J. As the first invariant 

that could distinguish spaces which are homotopic but not homeomorphic, it can be used to 
classify lens spaces [HUH [19]. Analytic torsion (or Ray-Singer torsion) is a smooth invariant of 
compact manifolds of odd dimensions defined by Ray and Singer |23[ [24] as an analytic counter- 
part of the Reidemeister torsion. They conjectured that the two torsions are equal for compact 
^ ■ manifolds. This Ray-Singer conjecture was proved independently by Cheeger [6] and Miiller 

[20J. Another proof of the Cheeger-Muller theorem with an extension was given by Bismut 
and Zhang [2] using Witten's deformation. In [T7j, the current authors introduced an analytic 
torsion for the de Rham complex twisted by a flux form. It shares many properties with the 
Ray-Singer torsion but has several novel features. In this paper, we review these developments 
and introduce a new twisted analytic torsion for the complex of invariant differential forms on 
a principal circle bundle. 

The paper is organized as follows. In Section 1, we review Reidemeister's combinatorial 
torsion, Ray and Singer's analytic torsion and the Cheeger-Muller theorem. In Section 2, we 
describe the analytic torsion of the de Rham complex twisted by a flux form, its invariance 
under the deformation of the metric and flux form when the dimension is odd, the relation 
with generalized geometry when the flux is a 3-form and the behavior under T-duality for 3- 
manifolds. In Section 3, we introduce a new twisted analytic torsion for the complex of invariant 
differential forms on the total space of a principal circle bundle twisted by an invariant flux 
form. We show that when the dimension is even, such a torsion is invariant under certain 
deformation of the metric and the flux form. Under T-duality, which exchanges the topology 
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of the bundle with the flux form and the radius of circular fibers with its inverse, the twisted 
torsions are inverse to each other for any dimensions. 

1. Reidemeister and Ray-Singer torsions 

In algebraic topology, several groups can be assigned to a topological space such as its 
fundamental group, homology and cohomology groups. For example, the dimension of p-th 
cohomology group H P (X, R) is, roughly speaking, the number of p-dimensional "holes" in the 
space X. These groups are invariants in the sense that if two topological spaces are the same, 
or more precisely homeomorphic, then the corresponding groups are isomorphic. So if we find 
two spaces X and Y with H P (X, R) ^ H P (Y, R) for some integer p, then we can conclude that 
X is not homeomorphic to Y. In this way, we can distinguish different topological spaces. 

However, these groups are invariant not just under homeomorphisms, but under continuous 
deformations, or homotopy. Two spaces can be homotopic but not homeomorphic to each other. 
As none of the above groups can distinguish two homotopic yet different topological spaces, 
we need additional invariants to achieve this. Reidemeister torsion is such an invariant that is 
invariant under homeomorphisms but not under homotopy. It is a secondary invariant in the 
sense that it is an element in a certain space constructed from the above groups (which are 
regarded as the primary invariants). 

To recall its definition, we consider a topological space of dimension n with a finite triangu- 
lation K. The simplicial structure gives rise to a cochain complex (C'(K),5); we will take the 
real coefficients unless otherwise indicated. Roughly, the Reidemeister torsion is the alternating 
product of the determinants of the coboundary operators 5 on various C P (K) (p = 0, 1, . . . , n). 
More precisely, we take the adjoint 5^ of 5 : C P (K) — > C P+1 {K) with respect to the inner product 
under which the p-simplices form an orthonormal basis. We define the Laplacians A p = 5' 5+55* 
on C P (K). Just like in Hodge theory, we have H P {K) = ker A p . If H P {K) = were true for 
any p, then all A p would be invertible, and we could define the Reidemeister torsion as the 
number [26j U d 

n 

r(K) = n(detA p )(-^ + V2. 

p=0 

To take into account of the non-trivial cohomology groups, we choose a unit volume element 
r] p of ker A p C C P (K) for each p. Thus, we have rj p £ det ker A p = det H P (K). We define the 
Reidemeister torsion 

n n 

r(K) = n(det'A p )(- 1 ) P+ V 2( g) ?? (-ir 

p=0 p=0 

as an element of the line det H* (K) = &) p= q det H P {K)(~~^ P constructed from the cohomology 
groups. Here, det' means taking the determinant of an operator in the subspace orthogonal to 
its kernel. 

We make two remarks here. First, in the above construction, the torsion is only defined up 
to a sign as the unit volume elements depend on the orientation on the cohomology groups. A 
more intrinsic way is to define torsion as a norm on the determinant line bundle [2]. Second, 
the cochains on K can take values in a local system that comes from an orthogonal or unitary 
representation of the fundamental group of K or equivalently, a flat vector bundle over the 
underlying topological space. The torsion defined would then depend on this data. 

It can be shown that t(K) is invariant under the subdivision of K and is invariant under 
homeomorphism. It also satisfy some functorial properties. Happily, it is not invariant under 
homotopy. A celebrated example is the lens space, which is the quotient of the 3-sphere by 
a finite Abelian group. In fact, the lens spaces are classified by their cohomology groups and 
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the Reidemeister torsion [91 [8j [19] . More recently, equivariant torsion has been used to classify 
isometries and quotients of certain symmetric spaces up to diffeomorphisms |14] . 

An alternative way to define the cohomology groups is by differential forms. If A is a smooth 
manifold, let d: Q P {X) — > i? p+1 (A) be the exterior differentiation of the space of p-forms, with 
d 2 = 0. The p-th de Rham cohomology is 

= ker(d: fP{X} ^ A P+1 P0) 
dR{ ' im{d: QP~ l {X) -)■ QP(X)) ' 

The de Rham theorem states that there is a natural isomorphism H^{X) = H P (K), where K 
is a triangulation of X. 

A natural question is how to represent the Reidemeister torsion analytically. Suppose X is 
compact, orientable and is equipped with a Riemannian metric. Then we define the Laplacians 
A p = d)d + dd) , where the adjoint is with respect to the inner product on fl p (X) given by the 
metric. If X is compact, then H^^X) = ker A p and is finite dimensional by Hodge theory. 
Taking the unit volume form rj p of ker A p C J2 P (A), the Ray-Singer analytic torsion is 

n n 

r R 5W = n(Det'A p )(- 1 ) P+ V 2( g) r? (-ir 
p=0 p=0 

as an element, defined up to a sign, of the line detH'^(X) [23]. The above construction 
generalizes easily if the differential forms are valued in a flat Hermitian vector bundle. By 
replacing d in the Laplacian with the flat connection, we obtain Ray-Singer torsion that depends 
on a flat vector bundle. 

The determinants of the Laplacians, and of many other operators on infinite dimensional 
spaces, are defined by regularization using zeta-functions. Let A be a self-adjoint semi-positive 
operator acting on a Hilbert space. Suppose the positive eigenvalues are listed as < Ai < 
A2 < • • • , taking into account of multiplicities. Since the eigenvalues Aj typically go to infinity 
as i — > 00, it does not make sense to consider their product in the usual sense. We define the 
zeta-function of A as 

00 

C4(s) = 5>- S = tr',4- 

i=l 

which is a function of one complex variable s. Here tr' means taking the trace in the subspace 
orthogonal to ker A. If A is elliptic, then by the heat kernel method, Ca(s) is analytic in s when 
5fts is sufficiently large and it can be extended meromorphically to the complex plane so that it 
is analytic at s = 0. We define 

Det'A = e-&< >. 

When A acts on a finite dimensional space, the above definition reduces to the finite product 
of the positive eigenvalues of A. 

The Ray-Singer torsion is invariant under the deformation of the Riemannian metric on X 
when dim A is odd and therefore is a smooth invariant of odd-dimensional compact manifolds. 
Furthermore, it satisfies the same functorial properties as the Reidemeister torsion. Ray and 
Singer [231 [24"] therefore conjectured that their torsion is equal to the Reidemeister torsion 
for compact odd-dimensional manifolds. The case for lens spaces was verified in |22| . The 
conjecture was proved independently by Cheeger [6] and Miiller |20j . when the local coefficients 
are given by an orthogonal or unitary representation of the fundamental group. The result is 
also true when the representation is unimodular [21] . In |2], Bismut and Zhang extended the 
Cheeger-Miiller theorem using the deformation of Witten |29| . They established the variation 
of the torsion under the deformations of the metric on the manifold of arbitrary dimensions 
and the Hermitian form on the vector bundle with an arbitrary flat connection. 
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2. Analytic torsion of twisted complexes 

Let X be a smooth manifold and H, a closed 3-form H called flux form, which has its origin 
in supergravity and string theory [27|, [Tj . We have an operator du = d + H A • acting on fi* {X), 
which squares to zero. So it can be used to define a twisted de Rham cohomology H'(X, H). 
The twisted de Rham complex is only Z2-graded and so is the twisted de Rham cohomology. Let 
fi°(X), Q 1 {X) denote the space of differential forms on X of even, odd degrees, respectively. 
Then the operator du acts between these two spaces. 

In fact, the above flux form can be a closed form of any odd degree with no 1-form component. 
In |17j . we introduced the analytic torsion of the operator du when X is compact, which we now 
assume. The main difficulty was that the twisted de Rham complex does not have a Z-grading. 
Formally, given a Riemannian metric on X, the twisted analytic torsion is 

r{X,H) = Det / ^ (x) (4d H ) 1 / 2 Det^ I(x) (4dH)" 1/2 r ? o® 

where rj^ is a unit volume element of H k (X,H) (k = 0,1). However, unlike the Laplacian 
Ah = d H dji + djjd H , the operator d H du on f2 k (X) is not elliptic and the usual heat kernel 
methods seem inadequate. Instead, we use some properties of pseudo-differential operators. Let 
Pj. be the projection in fi k {X) onto the image of d H . Since Pj, is a pseudo-differential operator, 
we have [28] 

MiVV) = ^-i + C0 + o ( s ), 

upon meromorphically continuing the left-hand side to s = 0. The coefficient c_i is essentially 
the non-commutative residue [301 E2] of P^. It turns out that since Pj. is also a projection, 
c_i = |30j. Consequently, the zeta-function of the restriction of d H du to the image of d H is 
regular at s = and its determinant is defined. It is also possible to include a flat Hermitian 
vector bundle in the definition of twisted analytic torsion. 

Just as the usual de Rham cohomology groups, the twisted counterpart H'(X,H) is also 
invariant under homotopies [T7]. The twisted torsion t(X, H) satisfies a similar set of functorial 
properties as in [23J. Moreover, when X is compact and dimX is odd, t(X,H) is invariant 
under the deformation of the Riemannian metric on X, the inner product or Hermitian structure 
on the flat bundle and the flux form within its cohomology class [17] . Therefore the analytic 
torsion t(X, H) is a secondary invariant in the same sense but in the twisted setting. 

When H is a 3-form on X, the deformation of the Riemannian metric g and that of the flux 
form H within its cohomology class can be interpreted as a deformation of generalized metrics 
on X [17] ■ Recall that a generalized metric on A" is a splitting TX © T*X = T+X © T_X such 
that the bilinear form 

{x + a,y + fi) = (a(y) + P(x))/2, 
where x,y are vectors fields and a,/3 are 1-forms on X, is positive definite (negative definite, 
respectively) on T±X and such that {T+X, T_X) = [11]. The subspace T+X is the graph of 
g + B, where g is a usual Riemannian metric and B is a 2-form. Given a generalized metric, 
there is an inner product on Q'{X) called the Born-Infeld metric [11] , It can be shown that the 
effect of deformation H (— » H — dB on torsion is equivalent to taking the adjoint of the operator 
du with respect to the Born-Infeld metric [T7j. This amounts to deforming a usual metric g to 
a generalized one. Thus analytic torsion should be defined for generalized metrics so that the 
deformations of g and of H are unified. 

We consider a special case when H is a top form on an odd-dimensional orientable compact 
manifold X. The cohomology class of H is [H] £ H top (X, R) = R. Then the twisted analytic 
torsion is [T7] 

r(X,H) = \[H}\T I{S (X) m ®r ] ^, 
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where trs(A") is the Ray-Singer torsion. This provides examples of twisted torsion for 3- 
manifolds, when H has to be a 3-form. The calculation of t(X, H) is based on the work of 
Kontsevich and Vishik [J5] on factorization of determinants in odd dimensions, although we 
hope that there are simpler methoda^ which can be useful in more general cases as well. The 
factor \[H\\ is also the torsion of the spectral sequence in [27J. It is not clear whether there is a 
simple relation between the twisted torsion and the classical Ray-Singer torsion in the general 
situation. 

Whereas Reidemeister's combinatorial torsion precedes Ray-Singer's analytic torsion, the 
simplicial counterpart of the twisted analytic torsion is still missing in the general case. This 
is because the cup product in the simplicial cochain complex is associative but not graded 
commutative. Let h be a simplicial cocycle that represents the same cohomology class (under 
the de Rham isomorphism) as H. Then (5 + h U • ) 2 = (h U h) U • may not be zero. However, 
the situation simplifies if the degree of H or h is greater than dim X/2, in which case h U h = 
for dimension reason. In particular, this condition is satisfied if H is a top-degree form. In fact, 
a twisted version of the Cheeger-Miiller theorem holds in this case [T7j . 

We consider the behavior of the twisted torsion under T-duality. Let T be the circle group and 
7r: X — > M, a principal T-bundle with Euler class e(X) € H 2 (M,Z,). Suppose H is a closed 
3-form on X with integral periods. By the Gysin sequence, there is a dual circle fibration 
7T : X — > M, whose fiber is the dual circle group T, with a flux 3-form H on X such that |4] 

ir*[H]=e(X), n*[H] = e(X) 

(modulo torsion elements). Thus T-duality for circle bundles is the exchange of background 
H-Rux on the one side and the Chern class on the other. We have the following duality on 
twisted cohomology groups [3]: 

H 5 (X,H) ^ H l (X,H), H l (X,H) ^ H 5 (X,H). 

Consequently, det H'(X,H) = (det H*(X, H))' 1 . When dimX = 3, since H is a top-degree 
form, we get [17] 

t(X,H) = t(X,H)- 1 

under the above identification. In general, it is not known whether there is such a concise 
relation. Instead, we will explore an invariant version of the twisted torsion and its behavior 
under T-duality in the next section. 

Our method applies to other Z2-graded complexes [18]. For example, suppose X is a complex 
manifold and H G Q°> l (X). Let d H = 5 + H A • . If dH = 0, then B 2 H = just as in 
the de Rham case. Using the same argument, we can define an analytic torsion of the twisted 
Dolbeault complex [18]. Alternatively, we can take H G Q^ 2 (X) with dH = 0. Let Q~™(X) = 
T(A p T 1 '°(X) A^T ' 1 ^)*). Then H A • : ft- p ' q (X) -> Q-^ p - l ^ q+2 (X) and again, B 2 H = 0. 
The cohomology of Bh contains information of the deformation of twisted generalized complex 
structures [16]. A torsion as a secondary invariant in this case can also be defined. 



For example, affirmative answers to the following questions on functions of one complex variable would 
provide a simpler calculation. Suppose a sequence of positive real numbers < Ai < A2 < ■ ■ • goes to infinity 
fast enough so that the zeta-function £(s) = IS absolutely convergent and hence analytic when 5Rs is 

sufficiently large. Suppose £(s) can be meromorphically continued so that it is regular at s = 0. We partition 
the set of positive integers to a disjoint union of finite sets Ik (k = 1,2,...). Let Ah = Ilig/ fc Then 
Z(s) — YlkLi S ls absolutely convergent and hence analytic when 5Rs is sufficiently large. Can Z(s) be 
meromorphically continued so that it is regular at s = 0? If so, is it true that Z'(0) = C'(0)? 
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3. T-DUALITY FOR CIRCLE BUNDLES AND ANALYTIC TORSION 



3.1. Analytic torsion for the complex of invariant forms. Consider a principal T-bundle 
7r: X — > M. Suppose X is compact and H is a T-invariant closed 3-form on X. We consider the 
^-graded complex (fi'(X) 1 - , dn) of T-invariant differential forms on X. At first sight, it seems 
difficult to define torsion, as the asymptotic expansion of the heat kernel with a group action [5] , 
and hence the pole structure of the corresponding zeta-function, is rather complicated. However, 
given a connection on X, the space f2'(X) T is isomorphic to Q*{M) © ft'(M). Through this 
isomorphism, operators on fi'(X) acts on sections of bundles over M. The torsion is then 
defined by the regularized determinants of elliptic operators on M. 

Suppose the metric gx on X is T-invariant and such that the length of every circular fiber 
is equal to some constant r > 0. Then 

gx = TT*g M + r 2 A A, 

where gM is a metric on M, A G Q 1 {X) is a connection 1-form with the normalization J T A = 1 
and stands for the symmetric tensor product. Since any uj G Q k (X)^ can be written uniquely 
as ui = r l l 2 n*uji+r~ l l 2 Af\it*oj2 with uj\ G Q k {M) and 0J2 G Q k ~ 1 (M), there is an isomorphism 

<t>: Q k {X) T -> Q k {M) © Q k ~^(M) 

for k = 0, 1 defined by <j){oS) = (wi,Ct>2). 

Lemma 3.1. The above isomorphism <fi is an isometry under the inner products on Q'{X)" S 
and Q'{M) defined by gx and gu, respectively. 

Proof. We have * x w = r~ x l 2 A A tt*(*m oji) + r 1 / 2 7r*(*Af 0J2) and thus 



CO A *x ^ = / 0JiA*m^1+ / W2 A*M^2. 
X JM JM 

The result follows. □ 

Let dfc be the restriction of the operator du to Q k {X)^ and let d^ = (j>dj,(j)~ l . If we write 
H = ir*H 3 + A A tt*H 2 with H 3 G /2 3 (M) and F 2 G tt 2 (M) and denote by F G ^2 2 (M) the 
curvature 2-form of A (that is, tt*F = dA), then 



dm r~ l F 



y rH 2 - d H3/ 

on fl k {M) © ,(2 fe_1 (M). Since <fi is an isometry, we have dt = 0<i^0 _1 and A^ := d^d^ + 

4=t4=t=/ a ^ _1 - 

Clearly, A^ is a second order elliptic operator on M whose leading symbol is the same as the 
(untwisted) Laplacian. The projection Pj, = cpP^cj)" 1 onto the closure of im(d|) is a pseudo- 
differential operator on M of degree zero. By the same argument as in Theorem 2.1 of [T7], the 
zeta-function 

is holomorphic at s = and hence we can define the determinant 

Det' T d\d- k := Det'44 = e"^ '^). 

The analytic torsion for the T-invariant part of the twisted de Rham complex is then defined, 
up to a sign, as 

vr(X,H,r) := (Det' T dld- ) 1/2 (Bet' T djdi) -172 ^ ® J/f 1 G det H'(X, H), 
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where the unit volume elements rj^ (k = 0,1) are as before, since ker(A^) = H k (X,H) are 
invariant under T. In fact, if % is the unit volume element of ker(A^), then since <f> is an 
isometry, % = (det(/>)(^). 

Theorem 3.2. If X is compact and dimX is even, then tj(X, H,r) is invariant under the 
deformations of gx in the class of T-invariant metrics such that the length of every fiber is r 
and the deformations of H in the space of T-invariant 3-forms representing the same cohomology 
class. 



Proof. If the metric gM is deformed along a path parametrized by u € R, then 

= ' -di = ~ [a ^ 

where a := ■ ^ ^he 3-form H is deformed in its cohomology class with a parameter v € R, 
let B £ Q 2 {X) be given by §f = -dB. Then 



dv 



where (3 = <f>fi<f> 1 and (3 = B A • . Since dimM is odd, the method of §3.1 and §3.2 of p2] 
shows that 

(Det; 4J ) 1/2 (Detv ^~ir 1/2 *7o ® 

is invariant under these deformations. In both cases, the isomorphism (j) remains fixed. So 
Tj(X, H,r) is also invariant. 

Finally, if the connection A is deformed with a parameter w 6 M, let C = — Then 

3 4 r~ 7 i 9 4 r~t Jti 

where 7 = (q r o C ')- Since dimM is odd, following the proof of Lemma 3.5 of [T7j, we get 

d_ 

dw 



^-loglDet'd^CDet'^)- 1 ] = 2 £ (-l) fc Tr( 7 Q s ), 



fc=0,l 

where Q^; = 4>Q^4>~ 1 is the projection onto ker(A^). On the other hand, the isomorphism 
(j) varies as w. If 4>^ is the isomorphism at w = 0, then ^[</> ° (0^) -1 ] = 7- Since r/^ = 
(det ^) _1 (%), following the proof of Lemma 3.7 of [17] . we get 

A (r?g ^-i) = _ £ (_!)* Tr( 7 Q s ) r/o ® r^ 1 . 

fc=0,l 

Therefore ti(X, H,r) is invariant under this deformation. □ 



3.2. Behavior under T-duality. Recall that if 7r: X — > M is a principal T-bundle and is 
a T-invariant closed 3-form with integral periods, the T-dual bundle tt: X — > M is a principal 
bundle whose structure group is the dual circle group T. It is topologically determined by 
ci(X) = n*[H]. We now explain T-duality at the level of differential forms. We have the 
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commutative diagram 

X x M X 



p / \ p 



X x 





M 

where I xj/I denotes the correspondence space. Then p*[H] = p*[H] G H 3 (X Xm X,Z). 
Choosing connection 1-forms A and A on the circle bundles X and X, respectively, the formula 

T(w) = [ e p * AA P* A p*u, uen*(X) 

gives linear map T: (l~ k {X) -> l? FrT (X), k = 0, 1. Similarly, we define 5: /2 5 pf) -> fl^X) 
by 

S(u) = I e- p * AA P* A p*uj, lu G (2 m (X). 
Jt 

We next explain the construction of the T-dual flux form H on X. Let tt*F = dA and 
tt*F = dA be the curvatures of the connections A and A, respectively. Since H — A A F is a 
basic differential form on X, we have 

# = AAtt*F-tt*Q 

for some I? G I2 3 (M). Define the T-dual flux H by 

if = 7T*F Ai-7T*/2. 

Then H is closed. Since 

d{p* A A p*i) = -p*H + p*H, 

we have 

T o dj{ = dfi oT, dn ° S = S o d^. 
Therefore T-duality induces isomorphisms on twisted cohomology groups 

U: H h (X,H) -> A; = 0,1 

with inverse S* [3], and there is an isomorphism 

detT*: det H m {X,H) 9* (det H'(X, H))^ 1 . 

We will relate the twisted analytic torsions under this identification. 

Given the Riemannian metric gx on X or a triple (gM,A, r), we define the T-dual metric on 



X as 



9 X =Tc*g M + r 2 AO A 



or given by the triple [gu , ^4> r 1 ) so that is T- invariant and the length of every fiber is r . 
We study the T-duality map on invariant differential forms. 

Lemma 3.3. Under the above choices of Riemannian metrics, 

T-.n h (X) T ^n^{X) f , S:Q ] '{X) f ^n i ^ I (X) T 

are isometries for k = 0, 1. 

Proof For any u = r^Vcdi + r^^A A tt*u 2 G ^(X) T , T(u) = r'^H*^ + r l l 2 A A 7r*wi. 
The result follows from applying formula in the proof of Lemma 13. II to both cj and T{oj). 5 is 
the inverse of T. □ 



Theorem 3.4 (T-duality and analytic torsion for circle bundles). In the above notations, we 
have, up to a sign, 

{detT*){T T {X,H,r)) = T f (X,Hy 1 € (det H'(X, H, r -1 )) -1 . 

Proof. We denote the restriction of to ]2 fc (X) T by dj.. Since T is an isometry, we have 
T o dj- = aj-— - o T and hence T o (d^d^) = {d^— j-f%+r) ° ^ follows that T isometrically maps 
the space of .ff-twisted even (odd) degree harmonic forms on X to the space of .ff-twisted odd 
(even) degree harmonic forms on X. So T maps the unit volume elements of H'(X,H) to 
those of H'(X, H) up to a sign. T also maps isometrically on other eigenspaces, preserving the 
(positive) eigenvalues. We deduce that 

Ct(s,4%) = Ct( s >4+i4+t) 
and the result follows. □ 

Acknowledgments. We thank the referees for helpful comments. 
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